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ABSTRACT 

The trajectories of dust particles ejected from a comet are affected by solar radiation pressure as a function of their ratios of radiation 
pressure cross section to mass. Therefore, a study on the orbital evolution of the particles caused by the radiation pressure reveals 
the physical properties of dust on the surface of the comet nucleus. In the course of NASA's Deep Impact mission, the ejecta plume 
evolved under the influence of the radiation pressure. From the evolution and shape of the plume, we have succeeded in obtaining 
« 0.4, where yS is the ratio of the radiation pressure to the solar gravity. Taking into account f3 « 0.4 as well as the observational 
constraints of a high color temperature and a small silicate-feature strength, dust particles ejected from the surface of comet 9P/Tempel 
1 are likely compact dust aggregates of sizes a 20yum (mass ~ 10"** g). This is comparable to the major dust on the surface of comet 
IP/Halley (~ lOyum) inferred from in-situ measurements and theoretical considerations. Since such dust aggregates with j8 ~ 0.4 must 
have survived on the surface against jets due to ice sublimation, the temperature of ice in the nucleus must be kept below 145 K, 
which is much lower than equilibrium temperature determined by solar irradiation and thermal emission. These facts indicate that 
9P/Tempel 1 has a dust mantle composed of 20yum-sized dust aggregates with low thermal conductivities ~ 1 erg cm"' K"' s"'. 
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1. Introduction 

Short-period comets are likely born in the outer solar system. 
They are originally composed of ices and dust particles. When 
the comets approach the sun, volatile components of their nuclei 
start to sublime. A jet of volatile gases pushes out dust having a 
large surface-to-volume ratio and selectively leaves heavy dust 
of a small surface-to-volume ratio (Prialnik et al., 2004, for a re- 
view). Therefore, relatively compact particles form the so-called 
dust mantle of comets. Dust particles emitted from comets give 
us important information on the surfaces of the comets related 
to their histories. Such a dust particle feels solar radiation pres- 
sure and then its orbit is determined by the ratio /3 of the radi- 
ation pressure to solar gravity, which depends on the morphol- 
ogy, size, and composition of the particle. The yS values of par- 
ticles are derived from the the sh a pes of comet tails arid trails 
(Finson & Probsteml Il968t iFuUel 1 19891: llshi^uro et all l2007t 
[ishi guro. 2008.) . Note that tail and trail particles of masses 
10"^-10"^g occupy the largest end of the size distribution of 
cometary comae, since smaller particles are blown away by 
the radiation pressure after leaving the comet nuclei. In-situ 
measurements of dust in the coma of comet IP/Halley showed 
that a typical coma dus t particle has a mass of 10"'°-10"''g 
jKolokolova et al.L l2007h . Because large dust particles in tails 
and trails are most likely a minor component on the cometary 
surface, we cannot derive the p values of typical dust particles 
on surfaces of comets from tail and trail observations. Therefore, 
the determination of the /? value from observations of comet co- 



mae is preferable. iHavward et al.l (l2000h attempted to derive the 
j3 value from barred structures formed in a coma of comet Hale- 
Bopp(C/1995 Ol) and found /? < 1. Owing to the fact that a 
continuous anisotropic ejection of dust particles in a coma blurs 
the trajectory of the particles, it is difficult to tightly constrain 
the /3 value. 

The advent of a promising occasion for a tight constraint 
on the y6 value was the Deep Impact (DI) mission, which suc- 
cessfully colhde d a 366 kg impactor with Comet 9P/Tempel 1 
(lA'HearnL l2007h . The resulting ejecta plume was observed by 
space telescopes and by telescopes on the Earth, which emerged 
after the DI collision, expanded in about a day, and then vanished 
in several days. The plume produced by such a single collision 
contained typical dust on the cometary surface and the dissipa- 
tion of the plume reflected the blow-away of the dust due to the 
radiation pressure. The /? ratio is estimated to be 0.08-1.9 based 
on data analyses of a temporal evolution on the brightness peak 
position of the plume at the direction to the sun or the shape 
of the plume at a ce rtain ti me (iMeech et al.L l2005t iMilani et al.L 
l2007l IWaikeretalL l2007t iBoehnhardt et al.L l2007h . These esti- 
mates of j8 unfortunately include large uncertainties, because the 
analyses were based on limited observational data. 

In this paper, we succeed for the first time in obtaining a tight 
constraint on the p value of typical dust particles on the surface 
of the comet, by taking into account both the temporal evolution 
in the distance to the plume leading edge at the solar direction 
and the plume shape at the time when the shape is significantly 
modified by the radiation pressure. We then use the obtained y6 
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values to discuss the physical properties of dust particles, which 
were originally on the surface of the nucleus and produced by 
the DI impact. The successful determination of the fi ratio en- 
ables us to estimate the mass and size of the dust particles based 
on the observational and theoretical considerations and the tem- 
perature of cometary ices, which has been kept low due to the 
existence of a dust mantle on the nucleus. We will show that the 
thermal conductivity and porosity of the dust mantle could also 
be deduced from our estimate of the /? value. 

2. Model 

Dust particles were ejected from the nucleus of comet Tempel 1 
by the DI event. These ejecta made an ejecta plume. According 
to an impact cratering theory, higher velocity ejacta come from 
a shallower layer of the target surfa ce and ejecta from the deeper 
region have lower velocities (e.g., ICrofi [l980h . The ejecta in 
the leading edge of the plume are thus composed of dust parti- 
cles on the surface of the comet nucleus, while the ejecta from 
the deeper region are mainly in the inner plume. Therefore, we 
could obtain the information of dust particles on the surface of 
cometary nucleus from the leading edge of the plume. 

The motion of the ejecta is mainly governed by radiation 
pressure from the sun. The observation from Earth showed that 
the ejecta plume expanded isotropically around the impact di- 
rectio n within several hours from the DI event dFeldman et al.L 
I2007h . We assume that particles ejected from the surface just af- 
ter the DI collision have a velocity projected to the plane perpen- 
dicular to the line of sight, vq. It should be noted that vp i s lower 
than t he real ejecta velocity by a factor 1 .5-2 (llpatov & A'Hearnl 
We treat their orbits on the plane. We set the x-axis in the 
direction to the sun on the plane and the y-axis perpendicular 
to the X-axis. The position {x,y) of an ejecta at time t after the 
collision of DI impactor is given by 

X - tVQ cos J - f-pCM^ sin a/2D^, (1) 
y = fVQsiny, (2) 

where D ^ 1 .5 AU is the distance of Comet 9P/Tempel 1 from 
the sun at the DI event, a ^ 48° is the angle between the line 
of sight and the incident radiation from the sun, G is the gravi- 
tational constant, and 7 is the angle between the direction of the 
projected ve locity and the x axis. The /? ratio is given by (e.g., 
iBums etal.L [T979) 



AttcGMqIu 

where Cp, is the radiation pressure cross section of the ejecta, m 
is its mass, Lq is the solar luminosity, M© is the solar mass, and 
c is the speed of light. 

Although the ejecta caused by the DI have various values 
of j6 and vq, not all of them are observable. Jorda et al. (200^ 
investigated the ejecta plume profile adopting distributions of /3 
and vq ™d found that the total cross section of dust particles was 
almost determined by single values for vq and /3. We therefore 
assume dust particles have a single value of /3 and eject with a 
velocity vq. In Fig. \T\ the dashed line given by Equations ([TJ 
and (|2]i with certain /3 and vq is hereafter called the leading edge 
of the plume at a certain time after the DI collision, where we 
assume the ejection is isotropic. The plume profile is determined 
by the /3 and vq values. Since the radiation pressure modifies 
the plume shape, the distance of the plume leading edge from 
Tempel 1 depends on the angle - arctan(y/x); the distance 
at around » is shorter than that at a large \Q\. The spatial 



distribution thus arrows us to constrain and vp. Nevertheless, 
since the observed plume shape includes errors, fi and vo could 
not be determined solely from the profile at a certain time. The 
time evolution of the leading distance at = 0° could give a 
further constraint to fix these values. 




Solar 
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Fig. 1. Schematic for the ejecta plume leading edge profile. The 
profile is modified by solar radiation (dashed lines). For refer- 
ence, the profile without the radiation is represented by the circle 
with dotted lines. 



3. Results 

At first, we focus on the spatial distribution of ejecta plume. 
Hubble Space Telescope observations within several hours af- 
ter the DI colhsion showed that the plume of ejecta are almost 
isotropic for the front side of the impact within Q = -145°-35° 
(Feldma n et al.L l2007l) : we therefore assume an isotropic distri- 
bution of 7. The shape o f the plume was then significantly modi - 
fied after t ^ 1 day (e.g.. lTbzzi et al.Ll2007l:IWaIker et aUl2007h . 
IWalker et al.l (l2007h presented the distance of the plume front 
from the nucleus at / = 25 hours after the impact. In Fig. |2] 
we compare the model calculations (solid and dashed curves) 
giv en by Equations (H) and ^ with the observational data (dots) 
by IWalker et al.l (l2007h . Large /? leads to the substantial vari- 
ation of the distance at the solar direction, while the variation 
of V() produces the difference of the distance independent of 6. 
We have found several possible combinations of and vq to ex- 
plain the data: e.g., vq = 210 ms"' and /? - 0.4 01 /3 - 0.65 and 
Vq = 215 m s^' as shown in Fig.|2] However, there are still uncer- 
tainties of the values of vq and /3 determined only from the spatial 
distribution. Therefore, we try to give a further constraint on /? 
and Vo, by taking into account the time variation of the plume. 

Fig. [3] shows the temporal evolution of the distance to the 
leading edge to the observed ejecta plume at = 0° (symbols) 
in comparison with the model (solid curve). The relation be- 
tween the distance and time is linear due to a constant velocity 
until about one day, but the ejecta were then significantly de- 
celerated by the radiation pressure. The evolution predicted by 
Equation ([T) with /3 = 0.4 and vq = 210 ms"' is in good agree- 
ment with the distance evolution obtained from the observations. 
The time at the peak distance, fpeak, is determined by dx/dt - 0, 
given by 

^P'=* /3GM0 sin a 
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Fig. 2. Distance of the ejecta plume leading edge from Tempel 
1 as a function of angle to solar direction on the view plane, 
6 = arctan(3;/ji:). Dots represent o bservational data at 25 hours af- 
ter the Deep Impact collision by Walker et al.l (l2007h (the bright- 
ness contour with 2ADU). Equations ([Ij and (|2]l gives the ejecta 
plume front taking y6 = 0.4 and vq = 210 ms"' (solid curve) and 
yS = 0.6 and vq = 215 m s"' (dashed curve). The direction of ini- 
tial plume emergence is indicated by DI and that of the sun is by 
0. 



3f^)f ,)days. 
\ /3 /\210ms-i/ ^ 



(4) 



The distance decreases by the radiation pressure at f > fpeak 
and the ejecta plume is blown away completely at f ^ 2fpeak- 
If we choose /3 - 0.65 and vo = 215ms"', fpeak ~ 2 days; 
the leading distances after 2 days are inconsistent with the ob- 
servational data. In addition, the maximum projected distance 
in the sun-ward directi on obtained from 6 days of observations 
by iMeech et al.l (l2005h was 3 x lO^'km, which is explained by 
ejecta with p » 0.4 because the peak distance is given by 
vofpeak/2 ~ 2.8 X 104(y3/0.4)-'(vo/210ms-')2km. Therefore, 
taking into account the time variation and the spatial distribu- 
tion, j8 ss 0.4 and vq ~ 210 m s"' are the likely solutions. 



4. Discussion 

Several authors derived the /? values within 0.08-1.9 from the 
shape of the ejecta plume at a certain t ime or the t empo- 
ral evolution of brightness peak position (Meech et al.L l2005t 



Milani et al.', '2007'; 'Walker et al.', '2007"; 'Boehnhardt et al.', "200^ 
Richardson et al., 2007; Kelleyetal., 2009). Schleicher et al. 
(120061) estimated the lower limit of /3 through a sophisticated 
simulation using a Monte Carlo jet model. On the other hand, 
we have applied a simple model described by Equations ([T) 
and (|2]l, whose transparent dependences of /3 and vq easily 
give a constraint on them, and have analysed both the lead- 
ing distance evolution of the plume at the direction to the sun 
and the shape of the plume at 25 hours after the DI event. 
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Observations 

A Meech et al. (2005) 

• Sugita et al. (2005) 

O Boehnhardt et al. (2007) 

■ Furusho et al. (2007) 

□ Kadono et al. (2007) 

A Lara et al. (2007) 

^ Walker et al. (2007) 

V Bauer et al. (2007) 
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Fig. 3. Time ev olution of t h e lea d ing edge d i stance 
from ob servation s (Meech et al., 2005; Sugita etal., 2005 
Boeh nhardt et al.L i2007; Furusho et al.. 2007; Kadono et al 



2OO7I iLara et 311120071 IWker et all 120071: iBauer et al.L 120071) . 



The solid curve is given by Equation ([T]l with /3 = 0.4 and 
v'o = 210 ms"'. The dashed curve is with /3 = 0.65 and vo = 
215ms-i. 



We have found that dust mainly determining the shape of 
the ejecta plume has /5 x 0.4. Equation ([3J gives Cpr/m ^ 
5200 cm^g"' foryS = 0.4. The value coiTesponds to 0.6 fim in 
radius for spherical particles having interna l density 2.5 gcin_-^ 
with Cpr given by geometrical cross section (iBurns et al.Lll979l) . 
However, dust particles mainly composing the ejecta plume are 
likely compact dust aggregates, which account for the obser- 
vational data about a silicate-fea ture strength and a color tem- 
perature (lYamamoto et al.L[2008h . According to lYamamoto et al] 
(20081), compact aggregates are assumed to have fractal dimen- 
sion D = 2.5. We apply th e description of fractal aggregates 
given in iMukai et al.l (Il992h . using monomers of an aggregate 
with radius 0.1 /zm and interior density 2.5 g cm"^. We approx- 
imate C, 



^ nr^, where A is the projected area of an ag- 
gregate and - Tg V5/3 with rg being the gyration radius of 
an aggregate. Dust aggregates with /? = 0.4 approximately have 
rc ~ 40jum (m ~ 10"** g ). On the other hand, based on in- 
situ measurements of IP/Halley, compact dust aggregates with 
X 10/im (m ~ 10~'°-10~^g) contribute most t o the cross- 
section area on the surface of comet Halley ( McDonnell et al.L 
IT987llKoiokoiova et al.L I2b07h . Both DI ejecta from the surface 
of Tempel 1 and typical dust aggregates of comet Halley have 
size of the order of lOjum, which may be determined by the 
temperature history of icy core of the nuclei. Therefore, we con- 
clude that the surface layers of short-period comets are mainly 
composed of such large aggregates. 

A dust particle with the projected area A and mass m on the 
comet nucleus with mass M and radius R is blown away due to 
volatile s ublimation if APy $ GMm/R^, where Py is the vapor 
pressure dPrialnik et all 120041) . Since A/m » AttcGMqPILq from 



3 



H. Kobayashi et al.: The Dust Mantle of Comet 9P/Tempel 1: Dynamical Constraints on Physical Properties 



Equation O) with Cpr ~ A, the condition to blown out for par- 
ticles with y6 is rewritten as Py ^ MLq/A-kcR^Mq/S. Since water 
ice is the main volatile material, Py is determined by water ice 
sublimation. As shown above, particles ejected by the DI event 
mainly have p » 0.4, which survive against ice sublimation on 
the nucleus surface. From the fact, we can estimate an upper 
limit on the value of vapor pressure, which indicates that the 
temperature of ice in the comet nucleus should have been kept 
lower than about 145 K. The ice temperature estimated above is 
lower than that on the nucleus surface which is calculated from 
the equilibrium between solar iiTadiation and its thermal emis- 
sion. Such a low temperature of ice therefore implies that the 
nucleus is covered with a dust layer that has a low thermal con- 
ductivity. 

Finally we here derive the thermal conductivity of the dust 
layer The thickness of the dust layer, L^, is much smaller than 
the radius of the nucleus. The heat flux in the dust layer is 
given by k^iT^ - Ti)/Ld, where is the thermal conductiv- 
ity of the dust layer, Ts is the temperature at the nucleus sur- 
face, and T, is the temperature at the boundary between the 
dust mantle and the ice core. The surface temperature obtained 
form the fitti ng of a nucleus spect rum varies from 272 + 7K 
to 336 + 7 K (iGroussin et all l2007h . The excavation de pth due 
to th e DI impact indicates Ld ~ 100-200 cm ( Yamamot o et al.L 
12008 ). The conductive flux received on the surface of the ice 
core should equal the ice sublimation energy on the ice core sur- 
face; the equilibrium with Ti - 145 K results in ~ l-l((2"s - 
ri)/190K)-i(i'd/200cm)ergcm-' K"' s"'. Dust particles in the 
dust mantle of Temple 1 that we have obtained from /3 have a 
porosity of 0.99 (or a filling factor <p of 0.01). The mean free 
path of photons in an aggregate is small enough that the radiative 
conductivity is negligible in the interest temperature range. From 
laboratory experiments, the thermal conductivity was derived 
to b e 10^-10"^ erg cm~' K~' s~' for dust lay ers with (p = 0.1- 
0.6 (iKrause et al.Ll201 lHOundlach & BlumLl2012h . The relation 
^d <P'' derived from theoretical analysis can ex plain experi- 
ments for aerogels using power- l aw index p ^ 1 .5 ( Fricke et al.L 
Il992t iHrubesh & PekalaL 11994 |Lu et al.', '1995"). The thermal 
conduction for dust layers measured by Krause et al. (2011) is 
also well described by a power-law formula with p = 1.5 (al- 
though /:> X! 1.8 is better). From the fitting formula, we estimate 
the thermal conductivity to be 1-10 erg cm"' K"' s"' for dust ag- 
gregates with (p = 0.01. Dust aggregates forming the dust layer, 
which can explain several facts obtained from the DI event, have 
a high porosity and produce such a low thermal conductivity of 
the dust layer 
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